Численный анализ: теоретические основы и корректность постановки задачи

Прежде чем использовать мощь численных методов, таких как метод Рунге-Кутты четвёртого порядка или формулы Адамса-Мултона, мы должны задать фундаментальный вопрос: Существует ли решение, и оно ли стабильно? Теоретические основы задач с начальными условиями (ЗНУ) дают математическое «зелёное свет», гарантируя, что наши дискретизации сходятся к реальному физическому результату, а не к численному шуму.

Основа: непрерывность Липшица

Чтобы контролировать распространение ошибок, нам нужна функция $f(t, y)$, которая не «скачет» слишком сильно. Это формализуется через условие Липшица.

Определение 5.1: Условие Липшица

Функция $f(t, y)$ удовлетворяет условию Липшица по переменной $y$ на множестве $D \subset \mathbb{R}^2$, если существует постоянная $L > 0$, такая что:

$$|f(t, y_1) - f(t, y_2)| \le L|y_1 - y_2|$$

для всех $(t, y_1), (t, y_2) \in D$. Эта константа $L$ — «предел скорости» изменения функции по вертикали.

Пример 1: Анализ констант Липшица

Рассмотрим $f(t, y) = t|y|$ на $D = \{(t, y) \mid 1 \le t \le 2, -3 \le y \le 4\}$. По теореме о среднем значении (или свойствам абсолютных величин):

$|f(t, y_1) - f(t, y_2)| = |t|y_1| - t|y_2|| = |t| \cdot ||y_1| - |y_2|| \le |t| \cdot |y_1 - y_2|$.

Поскольку максимальное значение $t$ в нашей области равно 2, константа Липшица равна $L=2$.

Геометрическая целостность области

Мы не можем решать задачу с начальными условиями в области, испещренной дырами. Нам требуется выпуклость.

Определение 5.2: Выпуклое множество

Множество $D$ выпукло, если для любых двух точек $(t_1, y_1)$ и $(t_2, y_2)$ отрезок, определяемый выражением:

$$((1 - \lambda)t_1 + \lambda t_2, (1 - \lambda)y_1 + \lambda y_2)$$

для $\lambda \in [0, 1]$, также содержится в $D$. Это гарантирует, что ни одна часть траектории решения не «выйдет» за пределы допустимой зоны вычислений.

Теорема существования и единственности

Когда эти условия выполняются, мы применяем теорему 5.4: Если $f$ непрерывна и удовлетворяет условию Липшица на выпуклом множестве $D$, то задача с начальными условиями $y' = f(t, y), y(a) = \alpha$ имеет единственное решение $y(t)$. Это оправдывает методы, такие как метод Эйлера ($w_{i+1} = w_i + h f(t_i, w_i)$), или более сложные процедуры предиктор-корректор:

$WC = w_{i-1} + \frac{h}{24}[9f(t_i, WP) + 19f(t_{i-1}, w_{i-1}) - 5f(t_{i-2}, w_{i-2}) + f(t_{i-3}, w_{i-3})]$.

🎯 Основной принцип: корректность постановки задачи

Задача является корректно поставленной если существует единственное решение и оно непрерывно зависит от начальных данных. Если константа Липшица $L$ чрезвычайно велика, задача становится «жесткой». В жестких уравнениях переходные части быстро затухают, но их производные (величина $c^n e^{-ct}$) не затухают, что требует применения алгоритма 5.8: трапециевидный с итерациями Ньютона для поддержания устойчивости.

ВОПРОС 1

Какое условие строго необходимо для теоремы 5.4, чтобы гарантировать единственное решение задачи с начальными условиями?

функция $f$ должна быть дифференцируема повсюду в $D$

функция $f$ должна быть непрерывной и удовлетворять условию Липшица на выпуклом множестве $D$

константа Липшица $L$ должна быть меньше 1

область $D$ должна быть невыпуклой

ВОПРОС 2

Для задачи с начальными условиями $y' = 1 + t \sin(ty)$ на $0 \leq t \leq 2$, какова лучшая оценка константы Липшица $L$?

$L = 1$

$L = 2$

$L = 4$

$L = 8$

ВОПРОС 3

В контексте жестких уравнений, почему простые явные методы часто не работают?

Переходная часть затухает слишком медленно

Величина $n$-й производной $c^n e^{-ct}$ не затухает так быстро, как сама функция, что вызывает нестабильность

Нарушается выпуклость

Константы Липшица всегда равны нулю в жестких задачах

ВОПРОС 4

Что происходит с общей погрешностью, если мы значительно уменьшим размер шага $h$ за пределы значения $h = \sqrt{2\delta/M}$?

Погрешность продолжает уменьшаться линейно

Общая погрешность увеличивается из-за доминирования ошибок округления

Решение становится более устойчивым

Константа Липшица уменьшается

ВОПРОС 5

Как определяется метод «середины»?

$w_{i+1} = w_i + h f(t_{i+1}, w_{i+1})$

Замена $T^{(2)}$ в методе Тейлора на $f(t + h/2, y + h/2 \cdot f(t, y))$

алгоритма 5.8: трапециевидный с итерациями Ньютона

Производная формулы обратной разности Ньютона

Кейс-стади: численная корректность постановки и динамика

МАТХ007: Расширенное задание

Инженер моделирует снаряд массой $m = 0.11$ кг, выпущенный вертикально вверх. Динамика описывается уравнением $mv' = -mg - kv|v|$, где $g = 9.8$ м/с², $k = 0.002$ кг/м. Перед решением необходимо проверить математическую основу.

Задание 1

Используя закон Кирхгофа как аналог устойчивости: если имеются данные тока [3.10, 3.12, 3.15, 3.18, 3.24] для $t=1.00$ до $1.04$, приближённо найдите напряжение $\mathcal{E}(t) = L \frac{di}{dt} + Ri$ при $t=1.02$ с помощью центральной разности (шаг $h=0.01$) при $L=2$, $R=10$.

Решение:
1. Вычислим $di/dt \approx (i(1.03) - i(1.01)) / (2h) = (3.18 - 3.12) / 0.02 = 3.0$.
2. $\mathcal{E}(1.02) = 2 \cdot 3.0 + 10 \cdot 3.15 = 6.0 + 31.5 = 37.5$ вольт.

Задание 2

Докажите, что задача с начальными условиями $y' = y/t - (y/t)^2, y(1)=1$ на $[1, 4]$ корректно поставлена, и найдите фактическое значение при $t=2$.

Решение:
1. $f(t, y) = y/t - (y/t)^2$. $\partial f/\partial y = 1/t - 2y/t^2$. На ограниченной области, где $t \ge 1$, эта производная ограничена, что удовлетворяет теореме 5.3.
2. Фактическое решение $y(t) = t/(1 + \ln t)$.
3. $y(2) = 2/(1 + \ln 2) \approx 2/1.6931 \approx 1.1812$.

Задание 3

Объясните, почему жесткая система, как $u_2' = -24u_1 - 51u_2 - 9 \cos t + \frac{1}{3} \sin t$, делает явные методы, такие как метод Эйлера ($w_{i+1} = w_i + h \cdot f_i$), непрактичными.

Решение:
В этой системе собственные значения Якобиана имеют большие отрицательные вещественные части. Переходное решение затухает как $e^{-ct}$ при большом $c$. Для устойчивости явного метода Эйлера требуется $h < 2/|\lambda_{max}|$. Здесь $\lambda \approx -51$, поэтому $h$ должно быть чрезвычайно малым (примерно 0.03), что делает вычисления для больших $t$ неэффективными. Мы должны использовать алгоритма 5.8: трапециевидный с итерациями Ньютона вместо этого.